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Abstract 

This work deals with the focusing Nonlinear Schrodinger Equation 
in one dimension with pure-power nonlinearity near cubic. We con¬ 
sider the spectrum of the linearized operator about the soliton solution. 
When the nonlinearity is exactly cubic, the linearized operator has res¬ 
onances at the edges of the essential spectrum. We establish the degen¬ 
erate bifurcation of these resonances to eigenvalues as the nonlinearity 
deviates from cubic. The leading-order expression for these eigenvalues 
is consistent with previous numerical computations. 

I 


1 Introduction 

The focusing, pure-power. Nonlinear Schrodinger Equation for 'ip{x,t) G C, 

X e M"", f e R, 


idtfj = —A-p — I'^l^ (NLSp) 

hnds applications in quantum mechanics, optics, and other areas, and has 
seen intensive mathematical study in recent years (eg. [22l [I5]). (NLSp) 

AMS Subject Classifications: 35P15, 35Q55. 

Key Words: nonlinear Schrodinger equation, linearized operator, edge bifurcation, Birman- 
Schwinger formulation, resolvent expansion, Lyapunov-Schmidt reduction 
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famonsly exhibits solitary waves (sometimes called solitons), solutions which 
maintain a fixed spatial profile, and which are observed to play a key role 
in the dynamics of general solutions. One naturally asks about the stability 
of these waves, which leads immediately to an investigation of the spectrum 
of the linearized operator governing the dynamics close to the solitary wave 
solution. Systematic spectral analysis of the linearized operator has a long 
history (eg. [26l [12], and for more recent studies P HI |2ll [25]). 

The principle motivation for the present work comes from |1] where reso¬ 
nance eigenvalues (with explicit resonance eigenfunctions) were observed to sit 
at the edges (or thresholds) of the spectrum for the ID linearized NTS problem 
with focusing cubic nonlinearity. Numerically, it was observed that the same 
problem with power nonlinearity close to p = 3 (on both sides) has a true 
eigenvalue close to the threshold. In this paper we establish analytically the 
observed qualitative behaviour. Stated roughly, our main result is: 

for p ^ 3, p 3, the linearization of the ID (NLSp) about its soliton has 
purely imaginary eigenvalues, bifurcating from resonances at the edges of the 
essential spectrum of linearized (NLS 3 ), whose distance from the thresholds is 
of order (p — 3)"^. 

The exact statement is given as Theorem [4] in Section H] and includes the 
precise leading order behaviour of the eigenvalues. 

The eigenvalues obtained here, being on the imaginary axis, correspond to 
stable behaviour at the linear level. A further motivation for obtaining detailed 
information about the spectra of linearized operators is that such information 
is a key ingredient in studying the asymptotic stability of solitary waves: see 
p 0 p [m [201 [m [B [7] for some results of this type. Such results typically 
assume the absence of threshold eigenvalues or resonances. The presence of 
a resonance is an exceptional case which complicates the stability analysis 
by retarding the time-decay of perturbations. Nevertheless, the asymptotic 
stability of solitons in the ID cubic focusing NLS was recently proved in |10j . 
The proof relies on integrable systems technology and so is only available for 
the cubic equation. The solitons are known to be stable in the (weaker) orbital 
sense for all p < 5 (the so-called mass subcritical range) while for p > 5 they 
are unstable [niET], but the question of asymptotic stability for p < 5 and 
p 7 ^ 3 seems to be open. The existence (and location) of eigenvalues on the 
imaginary axis, which is shown here, should play a role in any attempt on this 
problem. 

The generic bifurcation of resonances and eigenvalues from the edge of 
the essential spectrum was studied by [ 8 ] and [23] in three dimensions. Edge 
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bifurcations have also been studied in one dimensional systems using the Evans 
function in |T8] and [19]. We do not follow that route, but rather adopt 
the approach of [8l |23] (going back also to |T^ , and in turn to the classical 
work [IE]), using a Birman-Schwinger resolvent expansion, and 

Lyapunov-Schmidt reduction. 

Our work is distinct from 0122] due to the unique challenges of working 
in one dimension, in particular the strong singularity of the free resolvent at 
zero energy, which among other things necessitates a double Lyapunov-Schmidt 
reduction procedure. 

Moreover, our work is distinct from all of [TH] CHI El EB] in that we study 
the particular (and as it turns out non-generic) resonance and perturbation 
corresponding to the near-cubic pure-power NLS problem. Generically, a res¬ 
onance is associated with the birth or death of an eigenvalue, and such is 
the picture obtained in 0 |23l UHl |19]: an eigenvalue approaches the essential 
spectrum, becomes a resonance on the threshold and then disappears. In our 
setting, the eigenvalue approaches the essential spectrum, sits on the threshold 
as a resonance, then returns as an eigenvalue. The bifurcation is degenerate 
in the sense that the expansion of the eigenvalue begins at higher order, and 
the analysis we develop to locate this eigenvalue is thus considerably more 
delicate. 

The paper is organized as follows. The problem is set up in Section [21 
In Section [3] we collect some results that are necessary for the bifurcation 
analysis. Section 0] is devoted to the statement and proof of the main result. 
The positivity of a certain (explicit) coefficient, which is crucial to the proof, 
is verihed numerically; details of this computation are given in Section O 


2 Mathematical Setup 

We consider (NLSp) in one space dimension: 

idti! = -dli) - ( 2 . 1 ) 

Here -0 = i’ixit) : M x M —)■ C with 1 < p < cx3. The NLS fl2.ip admits 
solutions of the form 


'4){x,t) = Qp{x)e"^ ( 2 . 2 ) 

where Qp{x) > 0 satisfies 

-q; -q; + Qp = o. (2.3) 
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In one dimension the explicit solntions 

sech^ 

of fl2.3p for each p G (1, cxd) are classically known to be the nniqne solutions 
of fl2.3p up to spatial translation and phase rotation (see e.g. |3]). In what 
follows we study the linearized NLS problem. That is, linearize fl 2 .ip about 
the solitary wave solutions fl 2 . 2 p by considering solutions of the form 

t) = {Qp{x) + h{x, t)) e*h 

Then h solves, to leading order (i.e. neglecting terms nonlinear in h) 
idth = {-dl + l)h - QP-^h - (p - l)QP-^Re{h). 

We write the above as a matrix equation 

dth = JHh 

with 

r fRe{h)'\ T-^ — ~ 

\lm{h)) ^ [l 0 

H ■= ^ 0 

• ^ 0 -dl + 1- QP-^ 

The above JH is the linearized operator as it appears in [1]. We now consider 
the system rotated 

idth = iJHh 



and hnd U unitary so that, UiJHU* = where <73 is one of the Pauli 

matrices and with H self-adjoint: 


0-3 = 



U 



H 


f-dl + 1 0 \ 

VO -91 + y 


1 fp + l 

2 \p-l 


p- l\ 
p + ly 



H + V+\ 


In this way we are consistent with the formulation of [H [23]. We can also 
arrive at this system, idth = a^Hh, by letting h = [h hdj from the start. 
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Thus we are interested in the spectrum of 

Dp := (TaiL 

and so in what follows we consider the eigenvalue problem 


CpU = zu, 2 ; G C, M G L^(M, C^). 
That the essential spectrum of Cp is 

^ess(Dp) ( CX), 1 ] U [1, 00) 


(2.5) 


and 0 is an eigenvalue of Cp are standard facts [ 1 ]. 

When p = 3 we have the following resonance at the threshold z = 1 [1] 



( 2 . 6 ) 


in the sense that 


C 3 U 0 = Mo, Mo £ L°°, Mo ^ for g < CX). 


(2.7) 


Our main interest is how this resonance bifurcates when p 7 ^ 3 but |p — 3| is 
small. As is natural we seek an eigenvalue of fl2.5p in the following form 


z = 1 — a > 0 . 


( 2 . 8 ) 


We note that the spectrum of Cp for the soliton fl2.4p may only be located on 
the Real or Imaginary axes [1], and so any eigenvalues in the neighbourhood 
oi z = 1 must be real. There is also a resonance at z = —1 which we do not 
mention further; symmetry of the spectrum of Cp ensures the two resonances 
bifurcate in the same way. 

We now recast the problem in accordance with the Birman-Schwinger for¬ 


mulation (pp. 85 of IH]), as in |8l|23]. For fl2.8p . fl2.5p becomes 

{(t^H — 1 + a^)u = 


The constant-coefficient operator on the left is now invertible so we can write 


M = — 


{asH - 1 + =: 


Set 


w=\Vo\'2u, Vo:=V^P=^^ 
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and apply |Vo |2 to arrive at the problem 
w = -K^,pw, := 

with 

p(a) _ f 0 ^ 

“ 0 (-aj+ 2 -«=)-■ h 


(2.9) 


( 2 . 10 ) 


We now seek solntions {a,w) of (I2.9p which correspond to eigenvalues 1 — 
and eigenfunctions |Vo|“ 2 tc of fl2.5|) . The decay of the potential and hence 
|Vb| 2 now allows us to work in the space = L^(M, C^), whose standard inner 
product we denote by (•,•). 

The resolvent has integral kernel 


J_^-a\x-y\ 

R^^\x,y) = ' 2 “ 


0 


0 

1 „-V2-a'-‘lx-yl 


2\f2^c 


for a > 0. We expand as 


= —.R-i + Ro + aRi + a^R 


a 


R- 


( 2 . 11 ) 


These operators have the following integral kernels 


R-i{x,y)={^ U ) ,i?o(x,|/) = 


\x-y\ 


-V2\x-y\ j 1 Rlix^y') 


2 V 2 


\x-y\^ n 
4 ^ 

0 0 


and for a > 0 the remainder term R^ is continuous in a and uniformly bounded 
as an operator from a weighted space (with sufficiently strong polynomial 
weight) to its dual, and the entries of RR{x,y) grow at most quadratically in 
\x — y\. We also expand the potential in e where e := p — 3 

= Vo + eVi + e^V2 + e^Vr, e := p - 3 (2.12) 


and 


w = - 


w = 


2 1 
1 2 

1 /I 1 




2 VI 1 


Qi 


2 1 
1 2 


it/T _ 1 + l ^/3-l 

“2V^/3-l x/3 + 1 


Q 2 


Qs- 


Vi = - 
Vr = 


.lA 1 
2 A 1 

iA 1 

2 li 1 




Q 2 


2 1 
1 2 

2 1 
1 2 


Qi 


Qr 
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Here we have expanded 


Ql \x) = Ql{x) + eqi{x) + e‘^q 2 {x) + e^qn^x) 


and the computation gives 



By Taylor’s theorem, the remainder term qR{x) satishes an estimate of the 
form |gK(x)| <C{1 + |a:|^) sech^(a:/2) for some constant C which is uniform in 
X and £ G (—1,1). We will henceforth write 

Q for Qa and for Ka^p. 

3 Some Preliminaries 

We study (I2.9jl . that is: 


+ l)i(; = 0 . 


(3.1) 


Using the expansions fl2.1ip and fl2.12p for and we make the following 
expansion 


Kae — — (A"_io + eK-ii + e‘^K-12 + 
a 

+ A'oo + Kq2 + £^Kr2 

+ ctA'io + aeKR3 
+ a^KR4^ 


(3.2) 


where Kr 4 ^ is uniformly bounded and continuous in a > 0 and e in a neigh¬ 
bourhood of 0, as an operator on L^(R, C^). 

Before stating the main theorem we assemble some necessary facts about 
the above operators. 

Lemma 1. Each operator appearing in the expansion fl3.2p for Ka^e is a 
Hilbert-Schmidt (so in particular hounded and compact) operator from L^(R, 
to itself. 
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Proof. This is a straightforward consequence of the spatial decay of the weights 
which surround the resolvent. The facts that |||Vo|“ 2 || < C/ sech(a;), and that 
|Vo|^/^ decays like sech(x), while each of Vq, Vl, V 2 , decay like sech^(j/), and 
Vr decays at worst like sech^(3a;/4) (say if we restrict to |e| < |) imply easily 
that these operators all have square integrable integral kernels. □ 

We will also need the projections P and P which are dehned as follows: 
for f ^ Lf let 


Pf-- = 


{v,f)v 




1 

0 


as well as the complementary P := 1 — P. A direct computation shows that 
for any / G we have 


A_io/ = -^Pf. (3.3) 

Note that all operators in the expansion containing i?_i return outputs in the 
direction of v. 

Lemma 2. The operator P[Kqq + 1)P has a one dimensional kernel spanned 
by 

Wo := \Vo\^Puo 

as an operator from Ran(P) to Ran(P). 

Proof. First note that by fl2.7p 

VqUo = CT 3 M 0 - Huo, [VoUo]^ = [mo]i (3.4) 


from which it follows that 

Pwq = 0, i.e. Wo G Ran(P). 

Then a direct computation using (13.4p . the expansion (13.2p . the expression for 
Ro, and integration by parts, shows that 

(Poo + l)rco = 2n 

and so indeed P(Poo + ^)Pwo = 0 . 

Theorem 5.2 in mi shows that the kernel of the analogous scalar operator 
can be at most one dimensional. We will use this argument, adapted to the 
vector structure, to show that any two non-zero elements of the kernel must 
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be multiples of each other. Take w E with {w, v) = 0 and P{Kqq + l)t(; = 0. 
That is {Kqo + l)t(; = cv for some constant c. This means 


|V'o| 2 i?oho|Vo| + w = c\Vo \2 


2 M where u = 

fw \ 


V «2 / 


. We then obtain, after rearranging and 

V “2 / 

expanding 

c “ I /k 1^ “ ylQ^iy) (2mi(i/) + My)) dy 

/r^xp {-V2 \x - y\) Q^{y){ui{y) + 2 u 2 {y))dy 
We now rearrange the hrst component. Expand 


Ui 

U 2 


1 

'2 


k - yMi.y)Mii.y) + My))dy 
1 r 


2 

1 

2 


' —CXD 

poo 


X - y)Q‘^{y){2ui{y) + U 2 {y))dy 
{y - x)Q‘^{y){2ui{y) + U 2 {y))dy 


and rewrite the hrst term as 


f [ Q^{y)Mi{y)+ My))dy + l [ 

^ J — oo ^ J —oo 


X 


Q‘^{y){2ui{y) + U 2 {y))dy + h- 


yQ iy)i‘^My) + My))dy 

X / yQ^{y){‘2ui{y) + U 2 {y))dy 


where 

^■=1 [ yQ‘^iy)Miiy)+ My))dy 

Jr 

and where we used 2 Q^ui+Q^U2 = 0 since (w, v) = 0. So putting everything 
back together we see 


Ml 

U 2 


c + b +J^{x-y)Q^{y){2ui{y)+U2{y))dy 

-y\)Q iy)iMy)+ ‘^My))dy 


2\/2 


We claim that as a; —)■ cxd 


Ml 

M 2 




c + 6 
0 
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Observe 

poo 

/ {x-y)Q^{y){2ui{y)+U2{y))dy 


<-L 

0 


\y - x\Q‘^{y)\'2ui{y) + U 2 {y)\dy 

\y\Q‘^{y)\2ui{y) + U2{y)\dy 


as X —)■ oo. Here we have nsed the fact that w E implies Q\2ui + M2I ^ 
and that \y\Q E As well, in the second component 


^-V^\x-y\i ^2 


Q iy){uiiy) + 2u2{y))dy 


^g-v/2. 


e^yQ^{y){ui{y) + 2u2{y))dy 


+ e 


V 2 x 


’ —CO 

^co 




Q {y){My) + 2 u 2 {y))dy 


and 


3-\/2a: 


< e 


J —OO 
-'J2x I 

:,-V 2 x 


e^yQ‘^{y){ui{y) + 2u2{y))dy 
^ e^^Q‘^{y)\ui{y) + 2u2iy)\dy 


< e 

< Ce-^^ 

< Ce-^^ 


\ 1/2 

e‘^^yQ‘^{y)dy] 

> / 


Q^{.y)\ui{y) + ‘2u2{y)\‘^dy 


1/2 


< 


fX \ 1/2 

e^'^yQ‘^{y)dy\ 

— OO / 

rx \ 1/2 

e^^ve-^vdy ) 

— OO / 

^ 0 , x ^ cx) 


where we again nsed Q\ui + 2 m2 | € L^- Similarly, 


3\/2a; 


e ^yQ‘^{y){ui{y)+ 2 u 2 {y))dy 


-E 0 


as X —)■ cxo which addresses the claim. 

Next we claim that if c + 6 = 0 in (13.51) then m = 0. To address the claim we 
first note that if c + 6 = 0 then m = 0 for all x > X for some X, by estimates 
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similar to those just done. Finally, we appeal to ODE theory. Differentiating 
fl3.5p in X twice returns the system 

m" = —2Q‘^ui — Q‘^U2 (3.6) 

u'' - 2u2 = -Q^ui - 2Q^U2. (3.7) 

Any solution u to the above with u = 0 for all large enough x must be identi¬ 
cally zero. 

With the claim in hand we finish the argument. Given two non-zero ele¬ 
ments of the kernel, say u and u with limits as x ^ oo (written as above) c + h 
and c + h respectively, the combination 

* c + h _ 

u = u - -u 

c + b 

satisfies fl3.5p but with u*{x) —)■ 0 as x ^ oo, and so u* = 0. Therefore, u and 
u are linearly dependent, as required. □ 

Remark 3. Arguments similar to the estimates in the above LemmalB show 
that for a > 0 and w & Lf solving (12.9p the corresponding eigenfunction of 
fl2.5p u = \Vo\~ 2 w is in and so the eigenvalue z = \ — oA is in fact a true 
eigenvalue. 

Note that Kqq, and hence P{Kqq + 1)P, is self-adjoint. Indeed, a direct 
computation shows that Vq = — |Vo| and so 

/^ 00 = |G)|^i?oGo|I/o|-5 
= — I Vo h -Ro I ho 12 

= I Vo| ^ VoPol Vo| 2 

= (^oo)*. 

As we have seen above in Lemma [U, thanks to the decay of the potential, 
PKqqP is a compact operator. Therefore, the simple eigenvalue —1 of PKqqP 
is isolated and so 


{P{Kqq + 1)P) ^ : {x, Wo}^ ^ {v, Wo}^ (3.8) 

exists and is bounded. 

Finally, with the above facts assembled, we are now in a position to state 
the main theorem. 
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4 Bifurcation Analysis 

This section is devoted to the proof of the main resnlt: 

Theorem 4. There exists £o > 0 such that for —Eq < e < Eq with e ^ 0 the 
eigenvalue problem fl3.1l) has a solution {a, w) of the form 

W = Wo + EWi + e‘^W2 + w ^ ^ ^ 

2 ~ (4-1) 

Ct — S 0.2 ^ 


where 02 > 0, wq, w\, W 2 are known (given below), and |d| < C\e\^ and 
||ti}||L2 < C\e\^ for some C > 0. 

Remark 5. This theorem confirms the behaviour observed numerically in 
for p 7^ 3 but close to 3, the linearized operator JH (which is unitarily eguiv- 
alent to iCp) has true, purely imaginary eigenvalues in the gap between the 
branches of essential spectrum, which approach the thresholds as p —)■ 3. Note 
Remark\^ to see that u = \Vq\~' 2 w is a true eigenfunction of fl2.5l) . In 
addition, the eigenfunction approaches the resonance eigenfunction in some 
weighted If space. Furthermore, we have found that , the distance of the 
eigenvalues from the thresholds, is to leading order proportional to (p — 3)^. 
Finally, note that a = E^a 2 + 0{e^) with 02 > 0 gives a > 0 for both £ > 0 
and e < 0, ensuring the eigenvalues appear on both sides of p = 3. 

The qnantities in fl4.1l) are defined as follows: 


Wo : = 



Pwi := 
Pwi := 

Pw 2 : = 
Pw2 := 


02 : — 


-K_iiWo 

-{P{Koo + l)P)~" (^^PKooK.uWo + PKoiWo^ 

^{K_nWi + K_i2Wo + Oi2{Koo + l)wo) 

-{P{Koo + l)P)~' (^^PKooK.nWi + ^PKooK.uWo 

+ ^PKoo{Koo + l)tao + PKoiWi + PK02W0 + a 2 PKiow^ 

\(wo,KooK_iiWi) - \{wo,KooK_i2Wo) - {wo,KoiWi) - {wo,Ko2Wo) 

(wo, Kiowo) + 4(tao, iCoo(^oo + l)w^o) 
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Remark 6. A numerical computation shows 

a2 ^ 2.53/8 > 0. 

Since the positivity of a 2 is crucial to the main result, details of this compu¬ 
tation are described in Section\B 

Note that the functions on which P{Kqq + 1)P is being inverted in the 
expressions for Pwi and Pw 2 are orthogonal to both wq and v, and so these 
quantities are well-defined by fl3.8p . The identity 

1 

(wo, -A'ooA'-iiWo + KoiWq) = 0 

has been verihed analytically. It is because of this identity that the 0{e) term 
is absent in the expansion of a in fl4.ip . The fact that 

1 1 C (2 

0 = {wo, -iLooA'-iiWi -iLooA'_i2Wo + -^Koo{Koo+l)wo + KqiWi 

-h K 02 W 0 + a 2 KioWo) 

comes from our dehnition of Q!2- 

The above dehnitions, along with fl3.3p . imply the relationships below 


0 = K_,oWo (4.2) 

0 = K_iiWo + K_ioWi (4.3) 

0 = K_iqW2 + K_iiWi + K_i2Wo + 02(^00 + l)'?^o (4.4) 

0 = P(A'oo + l)w^i + PKqiWq (4.5) 

0 = P(A'oo + P)W2 + PKqiWi PK 02 W 0 a2PKioWo (4.6) 

which we will use in what follows. 


Using the expression for a in fl4.1l) . our expansion fl3.2p for now takes 
the form 

Kas = ~ (A"_io + eK-ii -\- e‘^K-12 + 
a 

+ Kqo + £Koi + Kq 2 + e^Kji2 
P {oi2e‘^ -f- a)KiQ P {oi2£‘^ + ofjeKji^ -|- (q;2£^ + 

=: — (K-iq P eK-ii P e‘^K-12 P e^Kjn) P Kqq P sKi P aK 2 
a 

where iLi is a bounded (uniformly in e) operator depending on e but not a, 
while K 2 is & bounded (uniformly in e and a) operator depending on both e 
and a. 
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Fnrther decomposing 

w = I3v -\-W^ (IF, v) = 0, 

we aim to show existence of a solution with the remainder terms a, (3 and W 
small. We do so via a double Lyapunov-Schmidt reduction. 

First substitute (14.ip to (13.ip and apply the projection P to obtain 

0 = P{Kc,,e + 

= P{Ka,£ + l)(wo + £UJl + £‘^W2 + (3 v + W) 

= P{Kqq + Ijwo + eP{K()o + Iju?! + sPKqiWq 

+ £‘^P{Koo + P KqiWi + P Kq2Wq + £^ CX 2 P KiqWq 

+ P{Kqo + l){^v + IF) + aPKioWo + P {eKi + dZa) {/3v + IF) 

+ £^P {Kr2Wo + Kq2Wi + K 01 W 2 + £Kq2W2 + eKR2Wl + £‘^Kr2W2) 

+ {(y2^‘^ F (x)PKiQ^EWi + £^W2) + (ci2^^ F di)£PK r^(wq + £Wi + £^W2^ 

+ {a2£^ F aYPKRi{wo + ewi + £^W2). 

(4.7) 

Making some cancellations coming from Lemma O fl4.5p and (14.6p leads to 
-P{Ko, F 1):P1F = 

l3PKooV F dPiFioWo F P {eKi F dPa) {/3v F IF) 

+ e^P {Kr2Wq + iFo2Wl F KqiW2 F £^ 02^2 F £Kr2Wi + £^Kr2W2) 

+ (cr2£^ F d')PKiQ^EWi + £^ 1 ^ 2 ) F (q;2£^^ F ^)ePKr^(wq + swi + e^W2) 
+ (q! 2£^^ F ci)^P Kr^{wq + £Wi + E^W2) 

=: P(fF; £,«,/?). 

According to (13.8F inversion of P{Kqq + 1)P on P requires the solvability 
condition 

PqP = 0 , Po := 11 ^ 112 (^ 0 , •)wo, Po := 1 - F), ( 4 . 8 ) 

which we solve together with the hxed point problem 

IF = (-P(Poo F 1)P) Po-F(lF; E.a, ^) =: Q{W- e, d, /3) (4.9) 

in order to solve (14.7p . 
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Write 

'.= P (^f3KQQV + oiKiqWq + {eKi + diK-^ (/?n + W) + E^fi + Eaf2 + 

where /i and /2 denote functions depending on (and bounded uniformly 
in) e but not a, while hi denotes an function depending on (and uniformly 
bounded in) both e: and a. 

Lemma 7. For any M > 0 there exists eg > 0 and R > 0 such that for 
all —Eq < e < Eq {but e 7^ 0) and for all a and (3 with |d| < M|e:p and 
\(3\ < M|e|^ there exists a unique solution W G L^n{n, tco}'*' of 114.9 1) satisfying 
\\W\\L2<R\Ef. 

Proof. We prove this by means of Banach Fixed Point Theorem. We must 
show that Q{W) maps the closed ball of radius R\e\^ into itself and that Q{W) 
is a contraction mapping. Taking W & orthogonal to v and tOo such that 
||hP||i,2 < R\e\^ and given M > 0 where |d| < M|e:p and \f3\ < we have, 

using the boundedness of (—P(A'oo + 1)^ Po, 

< CmiPKoov + P {eKi + dPa) + C\a\\\P {KioWo + e/s + ahi) 1^2 
+ C\\P {eKi + dPa) W\\l2 + |e|='C'||P/i|U2 

< CM\e\^ + CM\e\^ + C'|e|||W|U2 + C'|d|||W|U2 + C\e\^ 

< C\e\^ + CR\e\^ 

< R\Ef 

for some appropriately chosen R with |e| small enough. Here C is a positive, 
hnite constant whose value changes at each appearance. Next consider 

\\g{Wi)-g{W2)\\L^ 

< C\\P {eKi + cfK^) \\l^^l4Wi - W2IU2 

< C\e\\\PKi\\l2^l4Wi - W 2 IU 2 + C\a\\\PK2\\L2^L4Wi - W 2 IU 2 

< C'|e|||Wi - W2IU2 < «||hPi - W2IU2 

with 0 < K < 1 by taking |e| sufficiently small. Hence g{W) is a contraction, 
and we obtain the desired result. □ 

Lemma [7] provides W as a function of a and (3, which we may then substi- 
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tute into (14.81) to get 
0 = {wo, J^) 

= (]{wo, Koov) + d(wo, Kiowo) + el3{wo, K^v) + af3{wo, K 2 V) 

+ e^{wQ, fi) + £d(wo, / 2 ) + d^(wo, hi) + e{wQ, KiW) + d(wo, K2W) 

=: (3{wq, Kqqv) + d(M;o, KiqWo) + J^i (4.10) 

which is the hrst of two eqnations relating a and /5. 

The second eqnation is the complementary one to fl4.7p : snbstitute fl4.ip 
to fl3.ip bnt this time mnltiply by a and take projection P to see 

0 = aP{Ka^e + 

= K_iqWq + e{K_iiWQ + K_iqWi) 

+ {K_ioW 2 + K_iiwi + K_i 2 Wo) + £^02(7^00 + 

+ e^{K_iiW2 + K- 12 W 1 + KriWq + eK_i2W2 + sKriWi + e^KRiW2) 

+ l3K_iiiV + K_ioW + e{K_ii + eK_i2 + e‘^KRi){(3v + IT) 

+ d(iCoo T l)'M^o T ^^ 0 ' 2 P{Kqq + l)(rci + £102) + EotPi^Kim + l)(rci + £102) 
+ £^a 2 P{Koo + l){Pv + W) + aP{Koo + l)(/5n + W) 

+ aP^EKoi + Kq2 + £^Kr2 + CX,Kiq + CxeKr^ + O^Kri) 

X (wo + EWi + £^W2 + /3n + IT). 

(4.11) 

After nsing known information abont wq, wi, W 2 , a 2 coming from fl4.2p . fl4.3p . 
fl4.4p and noting that iC_iolT = —4P1T = 0 from fl3.3p we have 

0 = l3K_iov + d(iCoo + l)tt'o 

+ £^(A'_iiM; 2 + A'_i2tCl + KriWo + £K_i2W2 + eKriWi + £^^^ 1 ^ 2 ) 

+ £{K-ii + £^-12 + £‘^h^Ri){l3v + IT) 

+ £^oi 2 P{^Kqq + l)(rci + £’1112) + £oiP{^Kqq + l)(r( 7 i + £102') 

+ £^a 2 P{Koo + l)(/3n + IT) + dP(Poo + l){/3v + W) 

+ CxP^eKqi + e"^ Kq 2 + £^Kr 2 + Q-Kiq + (XEKr^ + (X^Kri) 

X [wq + £Wi + £‘^W2 + (3v + IT). 

Written more compactly, this is 

0 =/3K_iov + a{Koo + l)wo 

+ 5^/4 + £K^(^(3v + IT) + + CiKii(^f3v + IT) + a^h2 
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where is a bounded (uniformly in e) operator containing e but not a, while 
Ki is a bounded (uniformly in e and a) operator containing both e and a. 
Functions and /s depend on e (and are uniformly L^-bounded) but not a, 
while the function h 2 depends on both £ and a (and is uniformly L^-bounded). 
To make the relationship between a and (3 more explicit we take inner product 
with V 

0 = /3(n, K_iov) + d(n, (iLoo + f^) 

+ e{v, K^{(3v + W)) + oie{v, f^) + a{v, K^{l3v + W)) + a^{v, ^ 2 ) 

=: /?(n, K_iqv) + d(n, (-Ago + l)'?^o) + -^ 2 - (4-12) 


Now let 


and rewrite fld.lOp and fl4.12p in the following way 

Ar = ( (^o,AioWo) {wo,Koqv)\ (oi\ ^ fJ^i\ 

(Aoo + l)w,) (n, A.igu) ) J 

which we recast as a fixed point problem 

C = /i-'(';^d-ba.ftE)- (4.13) 

We have computed 

0 16 \ 

16 - 32 ) 

so in particular, A is invertible. We wish to show there is a solution {a,/3) of 
fl4.13p of the appropriate size. We establish this fact in the following Lemmas. 
Lemmas [ 8 ] and |9] are accessory to Lemma [TOl 

Lemma 8. The operators and functions K 2 , K 4 and hi, h 2 are continuous in 

d > 0. 



Proof. The operators and function in question are compositions of continuous 
functions of a. □ 

Lemma 9. The W given by Lemma\^is continuous in ( for sufficiently small 

kl • 
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Proof. Let (di,/9i) give rise to Wi and let {a 2 , ^ 2 ) give rise to IL2 via Lemma 
[71 Take |di — d2| < <^ and |/5i — ^21 < S. We show that ||hLi — lT2||i,2 < C5 for 
some constant C > 0. Observing K 2 depends on d, we see 


= II (P(i^oo + l)^ ' Po\\L^^L 4 J^{W,Xi;e) - J^{W 2 X 2 ;e)\\L 2 


< c 


(/3i — I32)Kqqv + {ai — a2)KiQWQ + e{Pi — P2)Kiv 
+ eKi{Wi - W 2 ) + aif^iK2{ai)v - a2/32K2{a2)v + diiC2(di)hLi 


— 0:2K2{o!2)W2 + s{pii — tt2)/2 + cfihiiai) — 02 ^ 1 ( 02 ) 

< CT + C'|£|||Wi-W2IU2 

+ ||Q;iiC2(di)(hLi — 1 ^ 2 ) + (diiL2(di) — d2iC2(d2)) hL2||L2 

< CT+ C'|£|||Wi - W2IIL2 


L2 


noting that |ai| < M|ep. Rearranging the above gives 


||W^i - H^2 ||l2 < C5 


for small enongh |e:|. □ 

Lemma 10. There exists £0 > 0 such that for all —Eq < £ < Eq ^ 0) the 
equation fl4.13p has a fixed point with |d|, \(3\ < M|ep for some M > 0. 

Proof. We prove this by means of the Bronwer Fixed Point Theorem. We show 
that F maps a closed sqnare into itself and that F is a continnous fnnction. 
Take |d|, |/3| < Mjep and and so by Lemma [7] we have ||hF||2,2 < |£|^i? for 
some i? > 0. Consider now 

I|24-^II l-^il 

< ll^“^ll ^4\P\\{'>^o,Kiv)\ + |d||/3||(M;o,7^2^^)| + |eH(wo,/i)| 

+ |e||d||(M;o,/2)| + \a4{wo,hi)\ + |£||(wq, FilF)| + |q!||(m;o, F 2lF)| 

< CM\£\^ + CM^\£f + + CM\£\^ + CM^\£\^ + CR\£\^ 
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and 


ll^-'ll 1-^21 

< ll^“i (^^\^\{v,f4)\ + \e\\{v,K3{/3v + W))\ + |d||e||(n,/ 5 )| 

+ |d|Kn, K/^{f3v + W)) \ + \a\^\{v,h2)\ 

< + CM\e\^ + CR\e\^ + CM\e\^ + CM^\ef + CMR\e\^ + CM^\e\^ 

< C\e\^ + CM\e\^ < M\e\^ 

for some choice of M > 0 and sufficiently small |£| > 0. Here C > 0 is a 
constant that is different at each instant. So F maps the closed square to 
itself. 

It is left to show that F is continuous. Given r; > 0 take \ai — q. 2 \ <6 and 
lA “ Al < Let (di, A) give rise to Wi and let (d 2 , A) give rise to W 2 via 
Lemma [3 We will also use Lemma [8] and Lemma [9l Now consider 


-7^1 (di, A) - -^i(d2, A)| 

= £(A - A) Ao, Kiv) + diA Ao, K2{ai)v) - d2A Ao, A'2(d2)n) 

+ £(di -d 2 )(wo,/ 2 ) + di(woAi(ai)) - d 2 AoAi(A)) 

+ e{wo, Ki{Wi - W 2 )) + di(wo, A'2(di)Wi) - d2(wo, A'2(d2)lE2) 

<G5 + G||A(di)-A(d2)||L2 

+ C'llfLi - hL2||i,2 + C'||A'2(di) - A'2(d2)||L2^2.2 
<C5 < 


V 


\\A-^V2 

for small enough 6. Similarly we can show 


\T,{auei) - T,{a,,P^)\ <CS< ' 

for 6 small enough. Putting everything together gives |A(Ci) — A(A)| < rj as 
required. Hence F is continuous. □ 

So hnally we have solved both (14.7p and (H.llh . and hence (13.Ih . and so 
have proved Theorem 01 
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5 Comments on the Computations 

Analytical and nnmerical compntations were nsed in the above to compute 
inner products such as the ones appearing in the definition of 02 (14. Ih . It 
was critical to establish that 02 > 0 since the expansion of the resolvent 
fl2.10p requires a > 0. Inner products containing wq and/or Pwi but not Wi 
can be written as an explicit single integral and then evaluated analytically or 
numerically with good accuracy. For example 

{wo, K 02 W 0 ) + ^(wo, A'ooA'_i 2 Wo) 

= - ^ / \x-y\- 3 Q^{x)) 

X {Q‘^{y)qi{y) - qi{y) + ^Q‘^{y)(i2{y) - ^2{y) - ^Q^{y))dydx 

^\^-^\{ 2 Q\x)- 3 Q\x)) 

X {Q‘^{y)qi{y) - qi{y) + ^Q^{y)q2{y) - 2 ^ 2 ( 2 /) - jQ^{y))dydx 

= - [ Q^{y){Q^{y)qi{y)-qi{y) + ^Q^{y)q2{y)-4:q2{y)-'^Q‘^{y))dy 

Jr ^ 

- / Q^{y){Q^{y)qi{y) - qi{y) + ^Q^{y)q2{y) - 2^2(2/) - yQ^(2/))d2/ 

Jr ^ 

^ - 2.9369 


2 V 2 


where 


C 2 — 77 / Q^qi — qi + 3Q^g2 ~ 4^2- 
^ Jr 

To reduce the double integral to a single integral we recall some facts about 
the integral kernels. Let 

h{y) = -l [ \x-y\{4Q^{x)-3Q\x))dx. 

^ Jr 

Then h solves the equation 


h" = -4Q2 + 3q4 

Notice that —4Q^ + 3Q"‘ = —2Q‘^U\ — Q‘^U 2 where u\ and U 2 are the components 
of the resonance Uq (12.61) . Observing the equation (13.6|) we see that h = Ui+c = 
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2 — + c for some constant c. We can directly compute h{0) = —2 to find 

c = —2 and so h = —Q^. A similar argument involving fl3.7p gives 


[ e y^{2Q^{x) - 3Q'^{x))dx = U 2 {y) 

2v 2 Jr 


-Q\y)- 


Computing inner products containing Pwi is harder. We have an explicit 
expression for Pwi but lack an explicit expression for Pwi. Therefore we 
approximate Pwi by numerically inverting P{Kqq + 1)P in 


P{Koq + 


-PPooA'-iitco + PKqiw^ 


=■ 9 - 


Note that {g^v) = {g,Wo) = 0. We represent P(iCoo + 1)-P as a matrix with 
respect to a basis The basis is formed by taking terms from the 

typical Fourier basis and projecting out the components of each function in 
the direction of v and Wq. Some basis functions were removed to ensure linear 
independence of the basis. Let Pwi = Then 

Ba = b 


where Bj k = {(pj, (Poo + l)0fc) and bj = {(pj,g). So we can solve for a by 
inverting the matrix B. Once we have an approximation for Pwi we can com¬ 
pute P(Poo + l)Pwi directly to observe agreement with the function g. With 
an approximation for Pwi in hand we can compute inner products containing 
Pwi in the same way as the previous inner product containing wq. In this way 
we establish that 02 > 0. 
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